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We consider orientational transitions of lamellar phases under the influence of a spatially nonuni- 
form electric field. The transition between parallel and perpendicular lamellar stackings with respect 
to the substrate is investigated as a function of the system parameters. The dielectrophoretic energy 
and the energy penalty for having dielectric interfaces perpendicular to the field's direction are iden- 
tified as linear and quadratic terms in a free energy expansion in the dielectric constant mismatch. 
We find that if the dielectric constant mismatch Ae is smaller than some critical value Ae c , parallel 
lamellar stacking will be realized, no matter how large the voltage difference between electrodes is. 
At Ae > Ae c , perpendicular stacking will appear if the voltage is high enough. Nonuniform fields 
remove the in-plane degeneracy present in the more common uniform fields. We therefore calculate 
the energy of grains of different orientations. The torque acting on the grains leads to the preference 
of only one orientation. The results have direct implications to block copolymer orientation and to 
surface patterning on the nanometer scale. 



Introduction. In recent years we have seen a large ef- 
fort directed toward finding ways to control the phase- 
behavior and orientation of self-assembled structures PD, 
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12^ . \2l\ have proved very useful. The use of electric fields 
is especially appealing, as the field strength scales favor- 
ably with the system size. Spatially uniform electric fields, 
however, pose a long-lasting problem since the orientation 
of the assembled phases is not unique - the symmetry of 
the field means that all grain rotations in the plane are 
energetically equivalent. 

Spatially varying fields remove this degeneracy, and 
thus can be quite useful in alignment of various 
mesophases. The early experiments of Russell et al. (28j 
have employed nonuniform fields, but ever since then all 
research have been on uniform fields. It seems that now, 
when such spatially uniform fields have been well under- 
stood and exploited possibly to their full potential, it is 
time to come back to spatially varying fields. In this ar- 
ticle we focus on the most simple periodic structure - 
the lamellar phase, which is found under the influence of 
an electric field emanating from a "razor-blade" electrode 
design (see Fig. 1). The lamellae are made up of two 
different materials, A and B, e.g., diblock copolymers. In 
this example, the two polymers A and B have different di- 
electric constants, eb and Eb- In the following we assume 
ion-free polymers; alternatively, for ion-containing poly- 
mers, application of a quasi-static field in the frequency 
~ 1 kHz renders the ions immobile but leaves the elec- 
trostatic equations unchanged p7| . In spatially uniform 
electric fields, the lowest-order contribution to the system 
electrostatic free energy is quadratic in the permittivity 
difference of the two constituents, Ae = ea — £b- As is 
explained in detail below, an inhomogeneity of the field 
gives rise to a dielectrophoretic force which is manifested 
in a linear term in Ae, and this has significance to the 
orientation selection and to phase-transitions [29] . 



(a) 



(b) L„ 



electrodes 

r V 

V 



(c) 



(d) L ± defect 




FIG. 1: Schematic illustration of the system, (a) Two thin "razor- 
blade" electrodes are laid down on the substrate. The voltage dif- 
ference between them is V. (b) Parallel stacking Lu - lamellae lie 
parallel to the substrate, (c) If the voltage is sufficiently high, elec- 
tric field can overcome interfacial interactions and prefer a perpen- 
dicular stacking L± (lamellae are parallel to the field lines), (d) A 
defect - an unfavorable perpendicular morphology where lamellae 
are perpendicular to the field lines. In subsequent calculations we 
took the distance between electrodes to be 1 /mi, and the lamellar 
period is 100 nm unless otherwise indicated. 



At this point it should also be pointed out that nonuni- 
form electric fields are in general neither interfacial nor 
purely bulk ones. In the razor-blade geometry, the field 
is high close to the electrodes' edge. However, sufficiently 
far from the electrodes the field behaves like E(r) = V/wr, 
where r is the distance from the middle of the gap. Thus, 
the integrated electrostatic contribution to the energy 
scales like 1/r. This energy indeed decays, but very 
slowly, and it has an important contribution even very 
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far from the electrodes. 

We assume that the lamellae are rigid enough so that 
the electric field does not bend them. In the example of 
block copolymers, this corresponds to the so-called strong- 
segregation regime, where xN 3> 1. Let us verify the 
validity of this assumption. The elastic bending energy 
per unit volume is written as F c i = ^K/R 2 , Where K 
is the bending modulus, and R is the local bending ra- 
dius (inverse curvature). For strongly stretched lamellae, 
K = D-fAB, where D is the lamellar period (D ~ 100 
nm) and 7,45 ~ 10 mN/m is the A/B interfacial tension. 
On the other hand, the electrostatic energy per unit vol- 
ume is F cs = ^eE 2 , where e is the dielectric constant 
and E the local field, which cannot exceed ^100 V//im 
because of dielectric breakdown. Let us take this max- 
imum value, in this case F cs — 10 5 J/m 3 . Therefore, 
F c \ — F cs if the lamellae are bent with a radius of curva- 
ture of 0.1 fim. The same estimate relates to the stresses 
(forces) of course. Electric fields cannot bend lamellae to a 
radius smaller than ~ 0.1 /im. In the razor-blade system, 
at larger distances from the electrodes the field is weaker, 
and therefore the lamellae should stay flat as well. Since 
the fields we consider are typically much weaker, we do 
not expect bent lamellae in this electrode arrangement. 
The above reasoning does not hold for weakly segregated 
lamellae: These lamellae have a much weaker modulus K 
and therefore significant bending can occur. 

A lamellar stack can therefore have the basic configu- 
rations: parallel or perpendicular to the substrate (Fig. 1 
(b) and (c)), denoted as L\\ and L±, respectively. Note 
that in the parallel stacking the first layer at the substrate 
is half as thick as the others. A third state exists which 
we denote the perpendicular defect. Here the lamellae 
normals are not parallel to the electrodes' edges. Fig. 1 
(d) represents the highest energy of such defects. Weakly 
segregated systems exhibit lower energy defects, e.g., T- 
junction or grain boundary. In some experiments with 
weakly segregated block copolymers on preferential sur- 
faces, few lamellar layers are adsorbed preferentially on 
the substrate (mixed morph olog y) due to the long-range 
effect of surface ordering [17ll30|. As is mentioned above, 
this system is out of the scope of the current work, and it 
will be dealt with in a subsequent publication. 

A peculiar feature of nonuniform fields is that the L» 
state can be favored over the L± one even in the absence 
of specific interfacial interactions with the substrate. In 
order to understand this, consider first the distribution of 
electric field squared for two semiinfinitc planar electrodes 
in the x-z plane, with a gap of 1 /im between them. This 
distribution is shown in Fig. 2 for a medium with spa- 
tially uniform dielectric constant e. Clearly E 2 is very 
high close to the surface and, in particular, close to the 
electrodes' edge at x = ±0.5 ^m. The field is small far 
from the substrate, and therefore interfacial instabilities 
are not expected [31| ; this is true even more so since above 
the electrodes' edges the field at y — > 00 is actually parallel 
to the substrate and also to the polymer/air interface. 

Let us now assume without loss of generality that 
Sa > £b- As is well-known in the field of dielectrophoretic 
forces [32| , a material with large value of e is drawn to re- 
gions with high fields, whereas small-e material is repelled. 



Since the electric field is largest near the electrodes' edges, 
an Lu state can form with the A-material touching the 
substrate. However, the work of Amundson et al. has 
shown that there is also a free energy penalty for having 
dielectric interfaces perpendicular to the field's direction, 
and this penalty is absent in the L± state. Clearly, the 
orientation selection depends on the magnitude of ea — £b- 
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FIG. 2: Plot of E 2 (x,y) in the x—y plane, for the case where the 
dielectric constant e is uniform, and the electrodes are at x > 0.5 
/an (V = I V) and x < -0.5 (im (V = -\ V). The largest field 
is at the electrodes' edge, x = ±0.5 fim. E 2 is scaled by 10 10 and 
given in (V/m) 2 . 

The electrostatic energy of the system is given by an 
integral over all space, 

F cs =-\J e(r)E 2 (r) d 3 r (1) 

The dielectric constant e(r) is a spatially varying quantity. 
In this study it is a periodic function. In the L11 state, for 
example, it is given by 

{e + iAe if nd < y < nd + ~d, n = 0,1,2... 

(2) 
£-§Ae if nd+ \ < y < (n + l)d, n = 0,1,2... 

where e = \{£a + £b) is the average dielectric constant, 
and the period is d. The above equation simply repre- 
sents a square-wave in the y-direction, where e alternates 
between ea and e b ■ The dielectric constant can be defined 
similarly for the other stackings. 

Theory and Results. Figure 3 shows F cs for the Ly and 

L± stackings at a fixed value of e = 6 and varying values 
of the dielectric constant mismatch. The electrostatic en- 
ergy is calculated numerically for a system with electrode 
gap of 1 /xm. F es (As) (dashed horizontal line) is constant 
for the L± case, because the electric field between the elec- 
trodes is independent of As and e. On the other hand, 
in the L\\ case (solid line), F cs (Ae) decreases first before 
it increases. The decrease is due to the dielectrophoretic 
term, linear in Ae, while the increase is due to the penalty 
associated with dielectric interfaces perpendicular to the 
field lines, scaling like (Ae) 2 . 

Let us make a short but very general mathematical di- 
gression which will clarify the last point. Denote Erj(r) 
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the electric field which corresponds to a system of uni- 
form dielectric constant and a given electrode design (not 
necessarily the one in Fig. 1). Eo is derived from a po- 
tential "00 ( r ) satisfying the proper boundary conditions on 
the electrodes: Eo = — VV>o- Suppose now that the dielec- 
tric constant changes from its average value by an amount 
£i(r): e(r) = e + £i(r). This change in permittivity leads 
to a change in field: E(r) = E (r) + Ei(r). We may now 
write the integrand of Eq. Jl} in the following way: 

/es = -\eV 2 = -\ml - ~ [e{E% + 2e"E • Ex] 

- i [2 £l E • Ei + eEl] -ieiE? (3) 

The first term on the right is the electrostatic energy of 
the system with uniform average e, while the other three 
terms are the deviations from it. The second and third 
terms (square brackets) are the dielectrophoretic and "di- 
electric interfaces" terms, scaling like e\ and ef, respec- 
tively. Finally, the last term scales like e\, and is small if 
£i -C £. For the case where this last term is dealt with 
the interested reader is referred to [20(. 



£i(r) = Ae • c(r). Thus, c(r) is a dimensionless function 
containing the spatial variation of E\, and whose spatial 
average vanishes: (c(r)) = 0. For the square- wave exam- 
ple of Eq. ([2J, c = ±|. It then directly follows that 

Ae ~ 

ipx = — tpi (r;c(r), geometry, V) (5) 

e 

where i^x obeys the equations 
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FIG. 4: Electrostatic energy F CB of perpendicular-defect structure 
(Fig. 1 (d)) as a function of Ae. F eB is normalized by its value when 
Ae = and is always increasing. Other parameters as in Fig. 3 



V 2 V'i = Vc • E 



(6) 



and ij) = on all electrodes. Since tpi is a universal 
potential independent of Ae, ipi is linear in £\je (and 
in fact it is linear in V as well). Similarly we find 
Ei = (Ae/e)Ei(r; c(r), geometry, V), with Ei indepen- 
dent of Ae. We now rewrite Eq. ([3]) as follows: 



FIG. 3: Numerically calculated electrostatic energy F CB [Eq. (I)] 
of parallel Ly (solid line) and perpendicular L± (horizontal dash- 
dotted line) stackings as a function of the permittivity difference: 
Ae 

— £ A & B • F cs is normalized by its value when Ae — 0. -Pes 
of perpendicular lamellae is constant, while that of parallel ones 
decreases before it increases (see inset). The critical value of Ae is 
Ae c — 1. When Ae < Ae c , Lu is preferred over L±. If Ae > Ae c , 
L± is preferred. We took the average dielectric constant to be e = 6, 
the lamellar period is 100 nm, and the electrode gap is 1 fim. The 
dashed line is a similar plot of F eB for Ln lamellae, with the same 
parameters; only the electrode gap is 2 fim. 

On the basis of this expansion and denoting Ei = 
— V-01, one can easily show that ip\ obeys the following 
equation 



V 2 0i 



1. 



:Vei • E 



(4) 



all 



with the boundary conditions that ipi = on 
conductors. Clearly ip\ can be written as ipi = 
0i (r, ei/e, geometry, V), where geometry refers to the 
electrode geometry and V to the electrode potential dif- 
ference (in the case of just two electrodes). We now 
write E\ in a form that puts emphasis on dimensions: 



(Ae) 2 

Ae/i + _ I 2 + const. 



(7) 
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c(r)E 2 + 2E • Ei 
2c(r)E • Ei + E 2 
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The expansion of -F cs is now transparent to order (Ae) 2 , 
as both Ji and I2 are independent of Ae, are quadratic in 
V 2 , and depend on geometry and c(r). In order to further 
demystify the above expansion, consider the simple one- 
dimensional example of uniform electric field Eq (parallel- 
plate capacitor), with c = ±i. In this case we find E\ = 
—cE and E\ = —c(Ae/e)E , and since (c) = we find a 
rather well-known result: (/ cs ) = |[(Ae) 2 /e]i? 2 + const. 

We now return to the razor-blade electrode design and 
the results presented in Fig. 3. The descent of i^s for 
parallel lamellae is due to a negative value of I\ , stemming 
from the dielectrophoretic force. The subsequent increase 
at larger value of Ae is due to a positive 12- The critical 
value of Ae, Ae c , is given by the relation 



Ae c 



-eh/h 



(8) 
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FIG. 5: (a) Illustration of a defect perpendicular morphology. 
Lamellae make an angle 8 in the x— z plane, as defined in (b). The 
system experiences torque which tends to align the stacking, prefer- 
ring the state with 8 = 0. 

The existence of Ae c is indeed important - at all Ae < 
Ae c the morphology is that of parallel layers irre- 
spective of the applied voltage or the magnitude of the elec- 
tric field. In uniform electric fields similar critical value 
of Ae does not exist. The value of the last term ignored 
in Eq. §3§ is numerically verified to be negligible in this 
calculation. 




e [deg] 

FIG. 6: Solid line: electrostatic energy F CB of perpendicular lamel- 
lae as a function of rotation angle 8 defined in Fig. 5. F es is 
scaled by |-F e s(f = 0)|. 8 = corresponds to "perfect" per- 
pendicular layering, while 8 = 90° is the defect with the high- 
est energy. The torque is L = dF CB /d8. Dashed line: a fit in- 
terpolating the maximum and minimum values by a sin 2 (8) fit: 
F CB = F CB (0) + [F cs (90°) - F es (0)] sin 2 (8). We took e A = 8 and 
eg = 4, yielding Ae = 4 and e = 6. The numerical accuracy for the 
point marked with a circle is questionable. 

In Fig. 4 we plot F cs as a function of Ae for the 
perpendicular-defect state sketched in Fig. 1 (d). At a 
given voltage and Ae, this state has the highest electro- 
static energy since the two electrostatic terms are unfa- 
vorable - the electrodes are not covered with the high-e 



material (ii > 0), and the field lines cross the lamellar 
interfaces {h > 0). 

Figure 5 depicts a lamellar grain in a defect state: the 
lamellae normals are not parallel to the electrodes' edges. 
The highest energy rotation has 9 — 90°, while the lowest 
is the state with 9 = 0. In Fig. 6 we present the elec- 
trostatic energy F os as a function of the rotation angle 9. 
The torque acting on the sample to orient it in the pre- 
ferred direction is given as the derivative: L = dF cs (9)/d9; 
it vanishes for the two extreme cases 9 = and 9 = 90° 
[l2l [13] • Indeed, when Eo is uniform in space we find 
F es (9) = F cs (0) + [F cs (90°) - F cs (0)} sin 2 (0). As is seen in 
the figure, the actual energy is higher than this estimate. 




V [Volt] 

FIG. 7: Phase diagram in the voltage-interfacial interactions plane. 
V is the voltage between the electrodes (see Fig. 1), and 7.4,3 and 
IBS are t ne interfacial interactions of the A and B polymers with 
the substrate. Above the solid line (green) and for Ae = 4, Ly is 
stable, while below it L± is expected. The dashed blue line is the 
same, but for Ae = 2. In both cases e = 4, Ae > Ae c — 1, d = 100 
nm, and the electrode gap is 1/im. 

Finally, the interfacial interaction of the two materials 
with the substrate can be taken into account as well. Let 
us call 7as and 7bs > 7as the interfacial energies per 
unit area of the A and B materials with the surface, re- 
spectively. The free energy difference between the L» and 
L± states is 

AF = I 1 Ae + I 2 ^^ + ^S( lAS - lBS ) (9) 

where S is the substrate area. The prevailing state is Lu 
if AF is negative and L± otherwise. On the basis of this 
free energy difference, one can construct a phase diagram, 
which is shown in Fig. 7 for two values of Ae. Note 
that both Ii and I 2 are proportional to V 2 , and since 
Ae > Ae c , the electric field terms favor the perpendicu- 
lar stacking. For fixed interfacial interactions, raising the 
voltage from small values to large ones destabilizes the 
.L|| and leads to perpendicular stacking L±. The critical 

voltage for this transition scales like (7AS — 7as) 1 ^ 2 - 

The polymer melt can be confined by another solid sur- 
face from the top. In this case there are two more 7^5 
and jbs corresponding to the second surface, and the aug- 
mented version of the equation above reads 

AF = J 1 Ae + 7 2 ^- + is(7ASi-7bBi) 
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+ t;S(jas 2 - 7bs 2 ) (10) 

where the "1" and "2" subscripts refer to the bottom and 
top surface respectively. Here we have assumed that the 
film is sufficiently thick so that the incommensurability 
between the lamellar thickness and the surface separation 
can be neglected and the parallel lamellae are not frus- 
trated, as is the case for surface separation larger than 
~ 10 lamellae. 

Conclusions. Lamellar phases under the influence of a 
spatially nonuniform electric field are considered. The 
role of the dielectric constant mismatch Ae is highlighted: 
the linear term in the free energy expansion is due to a 
dielectrophoretic force, while the quadratic term includes 
the free energy penalty for having dielectric interfaces per- 
pendicular to the field's direction. We have shown that a 
simple electrode realization which gives rise to nonuniform 
fields can bring about orientational transitions between 
several lamellar stackings. Specifically, for Ae < Ae c , par- 
allel lamellae are preferred over perpendicular ones even at 
very high voltages. When Ae > Ae c , there is an interplay 
between electrostatic forces and interfacial interactions. 



The "razor-blade" electrode design suggested here can 
find numerous applications in nanotechnology: the large 
torque is expected to remove the degeneracy between the 
L± states by orienting the lamellae perpendicular to the 
substrate and the electrodes' edges. More complex mor- 
phologies are expected to occur for block copolymers in 
the intermediate and weak segregations where the lamel- 
lar bending and grain boundary energies are smaller, and 
these systems should be systematically explored in this 
and more advanced electrode arrangements. 
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